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Background of the Invention 

1. Field of the Invention: 

This invention relates in general to radio frequency 
receivers subject to the reception of multipath .signals and in 
particular toreceivers measuring the time of arrival of signals 
received, such as those used in GPS navigation systems, which are 
adversely affected by multipath signals. 

2. Description of the Prior Art: 

GPS systems, as well as many other radio frequency (RF) 
communication systems utilizing frequencies high enough to be 
considered line of sight systems in which there must be a 
substantially direct line of sight between the transmitter (s) and 
receivers (s) for optimum operation, often suffer from multipath 



effects in which the receiver (s) must process signals received 
over a multiplicity of different paths. A common example is a 
simple broadcast TV system in which a TV receiver with an antenna 
receives multiple copies of the signal being transmitted. 

The multiplicity of signals being received results from 
additional, typically unwanted, signal paths including one or 
more reflections. When the signal path from the transmitter to 
receiver includes a reflection, this signal path must by 
definition be longer than the direct path. Multipath signals 
present a problem in systems, such as GPS systems, in which the 
time of arrival of the signal is to be measured or used because 
the time of arrival of the multipath signals depends on the 
length of the path(s) taken. 

The straightforward processing of all signals, including 
multipath or reflected signals, often degrades the processing 
performed by the receiver. The direct path is the shortest and 
therefore requires the least travel time from transmitter to 
receiver while the various unwanted multipath signals have 
various greater lengths, and therefore various longer travel 
times, than the direct path signals. 

GPS transmitters are positioned on satellites with complex 
orbital paths so that the locations of the multiple transmitters 
are constantly changing. This makes a highly directional antenna 
system almost completely unusable. Similarly, digital receivers, 
including those used in a GPS receiver, often do not rely solely 
on the amplitudes of the signals received, but rather rely on 
other signal characteristics, such as time of arrival. 

Multipath processing techniques currently used for complex 
receivers, such as GPS receivers, are often quite complex and 
subject to inaccuracies. An example of one such conventional 
technique is described in U.S. Pat. No. 5,414,729 issued on May 
9, 1995 to Patrick Fenton and assigned as issued to NovAtel 
Communications Ltd., Canada. In this technique, an 
autocorrelation function of a partially processed received 
signal, including multipath components, is compared to an 
estimated autocorrelation function of an estimated direct path 



signal to attempt to discern direct path signals from multipath 
signals for further processing. This technique of comparing 
processed and estimated correlation power, is complex and may be 
subject to error in that the partially processed signals relied 
5 on are themselves subject to degradation from many effects in 
addition to multipath effects including receiver limitations, 
which may reduce the accuracy or effectiveness of the multipath 
processing techniques. 

For example, in tracking a GPS C/A. signal to determine 
10 position information from GPS satellite transmitters, it is 

typically important to derive an accurate estimate of the time of 
j3 arrival, known as code phase, of the PRN modulation of the direct 
ifl path component of the C/A signals received from each of the 
^ various GPS satellites. It is also important to derive an 
10 accurate estimate of the phase of the underlying carrier signals 
^ transmitted from the satellites on which the modulation is 
; g applied, known as the carrier phase. However, as apparently 
shown for example in Fig.s 6, 7 and 8 of the above referenced 
Fenton patent, the delayed multipath components degrade the 
2JD tracking of the code and carrier phase estimates by distorting 

the correlation functions used is such tracking. 
Q What is needed is an improved technique for multipath signal 

processing which is less complex, less subject to error than 
conventional techniques and is applicable to a wide range of 
25 communication systems, signal encoding approach and conditions. 



Summary of the Invention 

30 In one aspect, the present invention provides an improved 

terrestrial navigation system using a GPS receiver in which the 
residual code phase tracking, or pseudorange, error due to 
simultaneous reception of multipath signals is detected, 
estimated and corrected. In particular, the distortion of the 

35 correlation function of the multi- and direct path signal 
composite as received with the internally generated code is 
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detected by comparison of an aspect of resultant correlation 
function with a model of the correlation function expected in the 
absence of multipath distortion. The comparison provides an 
indication of the sign of the residual error. 

It has been determined that the composite of a direct path 
and one or more multipath signals distorts the correlation 
function. If, as in most common cases, the multipath signal (s) 
are weaker than the direct path signals, the interference between 
such signals as received results in a predictable distortion of 
the correlation function. If the carrier phase (s) of the 
multipath signal (s) are shifted from between about 0° to 90° from 
the carrier phase of the direct path signal, the signals tend to 
reinforce each other resulting in a widening of the correlation 
function. Similarly, if the carrier phase (s) of the multipath 
signal (s) are shifted from between 90° to about 180° from the 
carrier phase of the direct path signal, the signals tend to" 
cancel each other resulting in a narrowing of the .correlation 
function. 

The correlation products are used in a code tracking loop to 
track and determine code phase. The most common scheme is to 
track points of equal magnitude (or power) separated by one C/A 
code chip width and estimate the time of arrival of the direct 
path signal as the mid-point between these points of equal 
magnitude. The points of equal magnitude on either side of the 
direct path arrival time are known as the early and late 
correlation time and the estimated arrival time of the direct 
path is called the punctual correlation time. In the presence of 
multipath signals, the correlation function has been found to be 
distorted so that the mid-point between the early and late 
correlations, that is the prompt correlation, * is not an accurate 
estimation of the arrival time of the direct path signal. 

When the correlation function of the composite is distorted 
to be wider than the correlation function expected for a direct 
path only signal, the distortion results in a lag erxor in which 
the prompt correlation lags the actual direct path signal 
received. Similarly, when the correlation function of the 



composite is narrower than expected, the distortion results in a 
lead error in which the prompt correlation leads the time of 
arrival of the direct path signal. 

The prompt correlation therefore leads or lags the actual 
time of arrival of the direct path signal by an amount designated 
herein as the residual code tracking error. The magnitude of the 
error can be approximated by the degree of narrowing or widening 
of the correlation function. Reduction of this error, or 
detection and correction of the error, enhances the accuracy of 
the resultant position determination. 

In another aspect, the present invention provides a spread 
spectrum receiver including means for correlating a received, 
code modulated spread spectrum signal with a code modulated 
signal replica at a first selected code phase delay and at a 
second selected code phase delay and comparison means for 
comparing characteristics of said correlations to derive 
information related to any residual code phase error in said 
second code phase delay. 

In another aspect, the present invention provides a method 
of processing coded, spread spectrum signals with potential 
multipath interference by generating a local replica of the code, 
correlating the local replica with coded, spread spectrum signals 
as received at a plurality of code phase delays and determining 
prompt code phase delay from differences in magnitude between at 
least two of said correlations. 

Brief Description of the Drawings 

Fig. 1 is an overview illustration of the operation of a car 
navigation system according to the present invention. 

Fig. 2 is a block diagram of the GPS car navigation system 
depicted in Fig. 1, used for improved navigation during reduced 
satellite visibility. 

Fig. 3 is a schematic representation of a single satellite 
channel of a GPS receiver used for fast satellite reacquisition . 



Fig, 4 is a schematic representation of a portion of the 
single satellite channel shown in Fig. 3 in which an additional 
plurality of sets of delayed code samples are correlated to 
provide a finer gradation of correlation intervals. 
5 Fig. 5 is a functional block diagram of a preferred 

implementation, on an ASIC, of the satellite tracking channels 
and associated processing components of the GPS car navigation 
system shown in Fig. 1. 

Fig. 6 is a functional block diagram of the Doppler Block of 
10 the GPS car navigation system shown in Fig. 1. 

Fig. 7 is a functional block diagram of the Coder Block of 
the GPS car navigation system shown in Fig. 1. 
ifl Fig. 8 is a functional block diagram of the Correlator Block 

/y s of the GPS car navigation system shown in Fig. 1. 

Fig. 9 is a function block diagram overview showing the 
My interconnections between the Doppler, Code, Correlator and other 
blocks of the system described in Fig. 5. 

Fig. 10 is a block diagram of the operation of the system, 
shown in Fig.s 5 and 9, illustrating the data path of the present 
2jfi invention. 

I;* Fig. 11 is a series of exploded time segments illustrating 

^ the operation of the data path of the present invention. 

Fig. 12 is a block diagram overview of a GPS receiver system 
illustrating a complete receiver system according to the present 
25 invention including a more detailed view of the satellite 
receiver section shown in Fig. 2. 

Fig. 13 is a graph of the correlation product of a direct 
path signal received without multipath interference, together 
with correlation products distorted by the presence of multi-path 
30 signals whose carrier phase differs from the carrier phase of the 
direct path signal by about 0° and by about 180°. 

Fig. 14 is a schematic block diagram of portions of a GPS 
receiver illustrating the delay-locked tracking loop with 
multipath residual code phase error detection, calculation and/or 
35 correction according to the present invention. 
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Detailed Description of the Preferred Embodiment (s) 

Fig. 1 is an overview illustration of the operation of a GPS 
5 car navigation system according to the present invention. The 
GPS car navigation system, described below in greater detail with 
respect to Fig. 2, is mounted in car 10 which is moving along the 
center of roadway 12. NAVSTAR satellite 14, in the lower left 
quadrant of the figure, is in view of car 10. A simulated GPS 
10 circular overhead display, positioned approximately over 

intersection 22 of roadway 12 and roadway 16 indicates that 
satellite 14 is between 0° and 45° degrees of elevation above the 
•Q horizon as viewed from car 10. 

For the purposes of illustration, satellite 18 is positioned 
lj| overhead between the elevation angles of 0° and 45° degrees. 
i|M However, the line of sight between satellite 18 and car" 10 is 
]q obscured by buildings 20 so that satellite 18 is not in view of 

car 10 at the position along roadway 12 as shown. Similarly, the 
;Jj line of sight between satellite 19 and car 10 is obscured by 
%Q buildings 21. However, as will be discussed below, when car 10 
jr crosses intersection 22, the line of sight between satellite 19 
12 and car 10, when the car is in position 11 within intersection 
22, may momentarily be clear. 

Turning now to Fig. 2, GPS car navigation system 24 is a 
25 first embodiment of a car navigation system according to the 
present invention which may be installed in car 10 of Fig. 1. 
GPS car navigation system 24 includes GPS car system module 2 6 
which is provided with signals received from satellites by GPS 
antenna 28, data related to .the then current - and expected 
30 future - physical environment of car 10 by for example Map Data 
Base 30 and data input from the operator of the car by for 
example input device 32. GPS car system module 26 provides 
output to the operator, for example, in the form of a GPS map 
display, via display unit 34, which may include both visual 
35 display as well as a voice interface announcing information as 
required to supplement or even partially replace visually 
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presented data. 

The present invention may be configured for use with only a 
GPS receiver, a GPS receiver aided by map data from, for example, 
Map Data Base 30 , and/or a GPS receiver aided by both a map data 
base as well as an external source of information, for example, 
from an external sensor. This external source of information may 
be used for maintaining position information by dead reckoning 
during those times when a sufficient number of satellites are not 
in view to provide the desired information. 

In operation, a composite of all signals received from 
NAVSTAR satellites is applied by GPS antenna 28 to satellite 
receiver section' 36 of GPS car system module 26. Signals from 
individual NAVSTAR satellites are then tracked in satellite 
specific tracking channels such as SatTRAK channels 38, 40, 42 
and 44. Although it is quite conventional to track 4 to 12 
satellites and therefore use 1 to 12 satellite tracking channels, 
only 4 such channels are shown herein for clarity. The outputs 
of these satellite specific tracking channels are processed by 
SatProcessor 46 to provide x use r, y U ser, z user and t use r data via 
appropriate logic control to a GPS position processor, such as 
PosProcessor or Nav Soln 48 which determines the navigation 
solution to determine position data. Position data is then 
applied by PosProcessor 48 to an appropriate display for the 
operator of the car, such as display unit 34. 

External sensor 49, in Fig. 2, may conveniently provide 
sensor data, or local or satellite position information, or 
position information which provided local position or satellite 
position information directly to PosProcessor 48 for comparison 
with the position information determined by SatProcessor 46 
and/or Map/Display Processor 50. External sensor 4 9 may 
conveniently be any sensor which provides information useful for 
updating position information for dead reckoning including 
direction, speed, velocity or acceleration or other data from 
which dead reckoning data may be derived. Conventional sensors 
include inertial navigation systems, with magnetic or optical 
gyroscopes, fluxgate compasses, odometer or wheel sensors or the 
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like. Alternatively, external GPS format signals, such as those 
provided by a pseudolite, may be used to update current satellite 
or position information. 

At the beginning of a navigated trip, the operator of car 10 
would typically provide data concerning the physical environment 
surrounding the intended route to GPS car system module 26 by 
inserting an appropriate data storage device such as a CD ROM, 
into Map Data Base 30 , and/or by entering the data via input 
device 32 which conveniently may be a keypad, keyboard, pointing 
device, track ball, touch screen, graphical pad, a voice 
recognition interface and/or a combination of such input devices. 

The operator of car 10 would also enter the intended destination 
into GPS car system module 2 6 via a data entry device such as a 
mouse or track ball interacting with display unit 34 and/or via 
input device 32. Map/Display Processor 50 of GPS car system 
module 2 6 would then develop the desired route, typically from 
the then current location as a point of origin to the desired 
destination, in accordance with the rules of navigation and 
details of the locale provided by Map Data Base 30. The 
appropriate route data is stored in Route Data Base 52, including 
the routing in the form of roadways and turns between roadways. 
Additional information, such as altitude, width of the roadways 
and etc. may also be contained within Map Data Base 30 and/or 
Route Data Base 52. These data bases may be contained within GPS 
car navigation system 24 and/or be made available to GPS car 
system module 26 from outside storage media such as diskettes 
positioned in appropriate disk drives. 

During navigation, each satellite in view may be tracked in 
a satellite tracking channel. If, for example, 4 or more 
satellites are in view, each of the satellites in view will be 
tracked in an individual channel, such as SatTRAK channels 38, 
40, 42 and 44. The output of the satellite tracking channels is 
then applied to SatProcessor 46 which would provide satellite 
based solutions of the four unknowns, such as x user , y user , z US er and 
tuser. The data represented by x user .and y use r are conventionally 
used as the two dimensional orthogonal components of the surface 
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of the earth such as north and east. However, in accordance with 
the present invention, x US er and y use r are preferably used to 
represent data for a pair of orthogonal directions specific to 
the direction of vehicle travel called the on-track and cross- 
track directions. 

Conventional bearing, such as north, south, east and west 
are relative to the magnetic or true north poles of the earth, 
while on-track and cross-track, as used in the present invention, 
are bearings made relative to the expected direction of travel of 
car 10 at any particular point in the route. For example, while 
a 90° turn from a heading of due north would change the angle of 
the vehicle velocity vector from 0° to 90° if bearings relative 
to the Earth's surface such as north and east are used, the same 
turn would show no change in the 0° angle of the vehicle velocity 
vector before or after the turn as long as car 10 remained on the 
expected track. 

The data represented by z user is typically surface elevation, 
such as the elevation above sea level, while the data represented 
by tuser is the exact time as determined from one or more of the 
satellite tracking channels. 

Solutions for all 4 unknowns of position information may be 
derived from signals from 4 satellites in view, so that exact 
position information within the limit of the accuracy then 
available from the GPS satellite constellation in view can 
therefore be applied by PosProcessor 48 to Map/Display Processor 
50. The position information determined from the satellites is 
processed with the physical data from Map Data Base 30, and/or 
the desired routing data from Route Data Base 52, to provide 
appropriate navigation information to the operator of car 10 via 
display unit 34 . 

If less than 4 satellites are in view, the t US er solution 
applied to PosProcessor 48 may be replaced by test 54 estimated 
solution derived for example from an internal clock model 54 in 
position estimate or model 63. Similarly, the Z US er solution may 
be replaced by z es t 56 solution derived from elevation estimate 
56, also in position model 63, in accordance with routing data 
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derived from Route Data Base 52 in accordance with then current 
GPS position information applied to Map/Display Processor 50. 
Test 54 and z es t 56 are applied to PosProcessor 48, and used in 
lieu of tuser and z user / when only two satellites are in view. The 
use of estimated or modeled solutions for the t and z variables, 
that is the use of t es t 54 and z es t 56 are conventionally known as 
clock hold and altitude hold, respectively. 

It must be noted that the particular configuration of GPS 
car navigation system 24 as described so far is only one of the 
many known ways of configuring such systems any of which may be 
used without departing from the spirit or scope of the present 
invention as defined by the claims. 

In accordance with the present invention, the width of the 
roadway, either known or estimated, may be used to provide y est 60 
for use in lieu of y US er when only one satellite is visible. y es t 
60 may be derived from Route Data Base 52 and/or Map Data Base 30 
Since the x and y unknowns are orthogonal, x user may be used to 
describe the on-track information, that is, the progress of car 
10 along its predetermined track while y es t 60 represents the 
cross track information, that is, how far car 10 has strayed from 
the center of the roadway. 

Referring therefore to Fig. 1, x use r is used to indicate the 
progress of car 10 along roadway 12 while y est 60 is used to 
represent the width of roadway 12. The actual width of the 
roadway may be derived from Map Data Base 30, or assumed because 
the actual value of the width of the roadway is relatively small 
and often therefore insignificant compared to the distances to be 
measured along the navigation route. Since the maximum allowable 
cross-track error, i.e. the maximum allowable appropriate value 
for y, is constrained by the physical width of the roadway, y est 
60 is relatively easy to accurately estimate. 

By using y es t 60, z est 56 and t est 54, it is possible to' 
provide useful navigation data for car 10 along a known roadway 
using signals from only a single satellite in view. It is 
important to note that reasonably accurate prior or initial 
position information may be required and that not all visible 
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NAVSTAR satellites will be suitable for single satellite 
navigation, depending upon the position of the satellite with 
respect to the path of car 10. The position information 
determined during single satellite navigation is along track 
position information which may be .accumulated and used for 
determining accumulated along track distance traveled. This data 
provides, and may be used in place of, the distance traveled 
information conventionally provided in a vehicle by an odometer. 

Referring now to both Fig.s 1 and 2, turning data may be 
used to improve terrestrial GPS navigation by using the detection 
of a known turn to update progress along a predetermined route. 
When at least 4 satellites are in view, the position of car 10 
may be known to the accuracy of the GPS system. When using 
clock, altitude or cross-track hold, or some combination thereof, 
the known position of the car is degraded by any inaccuracy of 
the estimate or estimates used. For example, during clock hold, 
internal clock model 54 drift and inaccuracy of the source of t es t 
54 will degrade the accuracy to which the position is known as a 
function of the magnitude of the inaccuracy. Similarly, any 
change in altitude from the estimated or fixed altitude, that is, 
any inaccuracy of z es t 56, will degrade the accuracy of the known 
position. Changes in roadway width and inaccuracies in the map 
data with regard to the roadway width, that is, any inaccuracy in 
Yost 60, may also degrade the position information. 

Even with 4 satellites in view, the geometry of the visible 
satellites may make it difficult to determine position by 
measurement of GPS signals. Further, during terrestrial 
navigation, it is not uncommon for satellites to be temporarily 
obscured from view during navigation by, for example, being 
blocked by buildings and other obstructions. 

It may therefor be desirable to update the accuracy with 
which the current position of the vehicle is known with actual 
position information whenever possible. The update information 
will sometimes be useful when 4 satellites are in view, but will 
always be useful as supplemental data when less than 4 satellites 
are in view. Update information is extremely useful during 
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single satellite navigation to avoid the accumulation of errors 
in position information. 

In operation, an original position and a destination were 
provided to the system which then determined the track to be 
5 followed. The track to be followed, or routing information, may 
be provided in the form of a data base of route information such 
as Route Data Base 52. In the example used, track 62 follows the 
centerline of roadway 12 to intersection 22 and then makes turn 
64 to follow the centerline of roadway 16. Track 62, roadways 12 

10 and 16, intersection 22 and turn 64 are provided to Route Data 
Base 52 during the preparation of the route by Map/Display 
Processor 50 from the then current position and the destination 

sfl entered via input device 32. 

t y s The physical position of car 10 is very accurately known 

In! when car 10 turns at turn 64. This accurate knowledge of the 
ip car's position at a particular time may conveniently be used to 
^ update the GPS navigation information by providing a position 

reset which is similar to a known initial position. Update 
!*! information from turns will most likely be useful if the angle of 
2B turn 64 is sufficiently large to provide an unambiguous position 
;™ determination. It is expected that any turn greater than 45° may 
1=5 be detected. As the speed of the vehicle increases, smaller turn 
angles may also provide useful information. The position update 
information is applied to position model 63 to update internal 
25 clock or t es t model 54, elevation or z es t model 56, y est model 60 as 
well as x est 61 which is a model of the along-track position of 
the car. These four estimates together for position model 63, 
which may be updated by information from Map Data Base 30, Route 
Data Base 52, current position processor 70, PosProcessor 48 
30 and/or external sensor 49, to form the most accurate available 

position model 63. Position model 63 may also be used to provide 
estimates to the same data sources. 

The actual turning of the car may be detected by a change in 
the vehicle velocity vector determined from the GPS data or from 
35 other conventional means such as a magnetic compass or an 
inertial navigation sensor. In accordance with the rapid 
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reacquisition system described below with respect to Fig. 3, GPS 
data alone may conveniently detect such turns even when single 
satellite navigation is required. The turn as detected by turn 
detector 66 is correlated with data from. Route Data Base 52 to 
5 determine the" actual position of the car to the accuracy of the 
Map Data Base 30. The accuracy of the data in Map Data Base 30 
may easily and conveniently be much greater than the accuracy 
available from the GPS system especially if single satellite 
navigation, or any combination of clock, altitude or cross-track 

10 hold, is used. Therefore the position update may provide a 

substantial increase in the accuracy of the then current position 

: ^ determination . 

:|g The benefit of the approach of this embodiment of the 

: ^ present invention is similar to the identification and use of a 
known waypoint during a dead reckoning navigation run. The 

Hp cumulative error is reduced substantially at the known waypoint 
so that additional, future position determination errors do not 
carry the burden of an accumulation of past errors. 

\Z As shown in Fig. 2, Route Data Base 52 provides data related 

j2k) to track 62, typically from Map Data Base 30, to Map/Display 
Processor 50 to display the current GPS position and may also 

;5 provide similar information to turn detector 66, turn comparator 
68 and/or current position processor 70 in order to update 
PosProcessor 48 with a position reset. 

25 Turn detector 66 may be configured in many different ways 

and is used to detect turns actually made by car 10 and select 
turns, such as turn 64, from Route Data Base 52 for later 
comparison with the detected turn. In accordance with a 
preferred embodiment of the present invention, turn detector 66 

30 may operate on the current GPS position provided by PosProcessor 
48 to develop a vehicle velocity vector position indicating both 
the direction and speed of travel. Substantial changes in the 
direction portion of the vehicle velocity vector would indicate a 
change in direction, such as a turn. Turn detector 66 may 
35 therefore detect turns directly from the GPS information by 

determining the vehicle velocity vector and detecting changes in 
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the vehicle velocity vector which represent a turn. 

Turn detector 66, or another unit if convenient, also 
operates on the route information provided by Route Data Base 52 
to determine the. expected position of car 10 along track 62 based 
5 on the then current GPS position information. Once the expected 
location of car 10 along the route is determined, one or more 
turns in the area of the expected position of car 10 can be 
selected for comparison with the indications of a physical turn 
derived from the GPS data, 
10 When changes in the actual vehicle velocity vector, as 

derived for example from the GPS position data, compare 
=2 appropriately with the changes predicted at a particular turn as 
;H derived from Route Data Base 52, the actual position of car 10 at 

the time of the turn can be very accurately determined and used 
Ip to update the GPS data at the turn. For example, if an actual 

turn is detected from a change in the vehicle velocity vector 
iQ from the GPS position of car 10 near the time predicted for that 

turn, the actual position of car 10 at the time of the turn can 
\ii be determined and used to update the then current GPS position 
24) for use as a position reset applied to PosProcessor 48. 
:Z Alternatively, turn detector 66 may use non GPS measurements 

=;3 for determining the occurrence of an actual turn of car 10, such 
as compass headings or inertial navigation determinations derived 
for example from external sensor 49, and applied directly to turn 
25 detector 66 or via PosProcessor 48 as shown in Fig. 2. • 
Detection of turns from GPS signals may easily be 
accomplished as long as 2 satellites are in view and provide 
appropriate geometries for determining two dimensional 
coordinates of the car's position. During single satellite 
30 navigation, as described above, the use of turn information for 
updating the last known position information becomes even more 
important, but the location of the single satellite in view, 
relative to track 62, becomes of even greater importance so that 
actual turns may be accurately detected. 
35 Turn detection may also be provided by monitoring changes 

between acquired and obscured satellites. If, for example, only . 
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satellite 14 was visible to car 10 on roadway 12 before 
intersection 22, and upon entering intersection 22, satellite 19 
suddenly became visible while satellite 14 was momentarily 
obscured, the change over from satellite 14 to satellite 19 could 
5 be used to indicate a turn in accordance with the data from each 
satellite. Using a rapid reacquisition scheme, as described 
herein below, the actual position at which the change of 
direction, that is, the position at which the switch between 
satellites occurs, can be sufficiently accurately determined to 

10 permit precise position update information at the turn. 

Similarly, turn comparator 68 may conveniently be 
^ implemented within another component of the system, such as 
iQ PosProcessor 48, Map/Display Processor 50 and/or SatProcessor 46, 
/? s so that a candidate turn may be selected from the route data for 
ij§> track 62 for comparison with the detected turn data. 
j V Referring now to Fig. 3, in another embodiment, the present 

invention provides for fast reacquisition of satellite signals, 

11 useful for example when a previously acquired satellite is 
obscured and then appears perhaps for only a short time, for 
example, as a car travels through an intersection. 

! ~ Referring to the line of sight between car 10 and satellite 

^ 19 as shown in Fig. 1, it is common in an urban environment for 
the buildings along the sides of. the street to act as a barrier 
wall obscuring the lines of sight to many GPS satellites. 

25 However, the barrier wall formed by buildings 20 and 21 is 

commonly breached at intersections such as intersection 22. For 
example, car 10 while traversing intersection 22 may reach 
position 11 in which the previously obscured line of sight to a 
satellite, such as satellite 19, is momentarily not obscured 

30 because of the break between buildings 20 and 21 at intersection 
22. This momentary visibility of a previously obscured satellite 
may occur while car 10 is in the intersection or at the edges of 
the intersection. 

The length of the momentary contact with satellite 19 is 

35 relatively short. For example, if intersection 22 is 60 feet 

wide and car 10 is traveling at 30 mph, the time taken to cross 
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the intersection may be as short as 1.3 seconds. Conventional 
GPS navigation systems would not reacquire and derive useful data 
from satellite 19, even if previously acquired, during this short 
time interval . 

5 In accordance with another embodiment, the present invention 

makes maximum use of such reacquisition opportunities by 
minimizing the time required for reacquisition, the collection of 
data and processing of the collected data for position 
determination. Referring now to Fig. 3, a portion of SatTRAK 
10 channel 38 is shown in greater detail as an example of the 

configuration of each of the satellite tracking channels. After 
original acquisition, SatTRAK channel 38 tracks a single 
! J3 satellite by operating on satellite signals 72 received by GPS 
; /: antenna 28. Satellite signals 72 include the signals from the 
ljj> satellite being tracked by SatTRAK channel 38 and are demodulated 
j |f and selected by being multiplied, in one of the correlators 74, 
= h by a copy of the 1023 chip pseudorandom, spread spectrum code 

applied to satellite signals 72 by the GPS satellite. 
;}rj Correlators 74 may be configured from exclusive NOR gates to 
W minimize the time required for providing a correlation result. 

During tracking, the copy of the code produced by code 
1=5 generator 76 and applied to exclusive OR correlators 74 by delay 
78 is synchronized with the code in satellite signals 72, as 
received, so that the copy of the code correlates with satellite 
25 signals 72. This may be accomplished in several different 
manners known in the art, including by shifting the time of 
generation of the code in code generator 76 and/or adjusting the 
amount of delay applied by an external delay. In any event, the 
code applied to exclusive OR correlators 74 when SatTRAK channel 
30 38 is locked to the selected satellite, is synchronized with the 
code being received from that selected satellite. This 
correlation is commonly called the on-time or prompt correlation 
to indicate this synchronization. 

Conventional GPS receivers maintain a lock on a satellite 
35 signal after acquisition by performing additional correlations, 
often called early and late correlations or correlations 
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performed by early and late correlators. These correlations are 
displaced in time by a certain delay such as one half the width 
of a C/A code chip from the on-time or prompt correlator. That 
is, if the time of occurrence of a particular chip in the 
satellite signals is time tO, the prompt correlator under ideal 
conditions would multiply satellite signals 72 with a replica of 
the code with the same chip at time tO. The early correlation 
would be performed at time tO - H chip and the late correlation 
would be performed at a time equal to tO + H chip. Whenever the 
synchronization between code generator 76 and satellite signals 
72 as received begins to drift, the correlation results begin to 
change in favor of either the early or late ' correlation at the 
expense of the prompt correlation. 

One conventional approach to maintaining lock on the signals 
from a particular satellite is to adjust the timing of code 
generator 76 with' a feed back loop used to maintain the power in 
the correlation products in the early and late correlators to be 
equal. In this way, code generator 7 6 may be continuously 
resynchronized with satellite signals 72 so that the accuracy of 
the system is within one half chip in either direction (early or 
late) of the signals received. 

When satellite signals 72 are temporarily lost, for example, 
because the satellite signals are temporarily obscured by 
buildings 20 and 21 as shown in Fig. 1, various techniques are 
used to attempt to synchronize code generator 76 with satellite 
signals 72 as received so that SatTRAK channel 38 can reacquire 
the signals from the desired satellite. As noted above, 
conventional techniques include clock and altitude hold and one 
embodiment of the present invention provides another technique 
called cross-track hold. 

However, unless the obscuration of the satellite signals is very 
brief, the accuracy of prediction of such techniques is not 
enough to maintain synchronization except for a very brief period 
of obscuration. 

In accordance with another embodiment of the present 
invention, massively parallel correlation is used to create an 
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expanded capture window of correlation capture around the then 
current predicted synchronization time in order to immediately 
reacquire a previously acquired, and then obscured, satellite 
signal. In particular, the speed of reacquisition is made 
5 sufficiently fast according to the present invention so that 
useful GPS position data may be acquired during the time car 10 
travels through intersection 22 even though, for example, the 
signals from satellite 19 were obscured by buildings 20 until car 
10 was within intersection 22. 
10 To this end, an expanded series of correlations are 

performed with a series of delays a fixed fraction of a chip 
width, such as H chip width, apart extending both early and late 
;|3 of the predicted prompt correlation. As shown in Fig. 3, 
*r s satellite signals 72 are devolved into a fixed number of samples, 
j]| by for example analog to digital conversion in A/D Converter 73, 
HU to provide n Signal Samples 75. A similar number of code samples 
are provided through k fixed H chip width delays 7 8 to provide k- 
1 sets of n Code Samples 80, progressing from a first set of n 
\Z Code Samples 80 with no delay to the k-lst set of n Code Samples 
2=8 80 which have been delayed by a total of k delays 78. It is 
j J convenient to use H chip delays for each delay 78, but other 
jSj fractions of a chip width may be used. 

The k/2th set of n Code Samples 80, or the set nearest k/2, 
may conveniently be delayed the correct amount to perform the 
25 prompt correlation in one of the exclusive OR correlators 74 with 
n Signal Samples 75 from A/D Converter 73 during tracking. The 
k/2th -1 set of n Code Samples 80 may then be used to perform the 
early correlation while the k/2th +1 set of n Code Samples 80 may 
then be used to perform the late correlation while tracking. The 
30 additional correlations may also be performed during tracking, 
but provide a substantial advantage when used during 
reacquisition . 

That is, in the present invention, the early, prompt and 
late correlations conventionally used in tracking may also be 
35 used during reacquisition mode, aided by a substantial number of 
correlations using additional delays. Whether or not the early 
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and late correlations are used, a convenient number of additional 
delays on each side of the prompt delay results from (k-l)=20 so 
that nine or ten H chip delays are provided on each side of the 
k/2th prompt delay. In this way, correlations are performed 
during reacquisition at time delays of 5 chip widths on either 
side of the predicted prompt or on-time delay. This represents 
an expanded capture window of on the order of ± 5 x 300 meters of 
potential error. That is, if the predicted synchronization with 
satellite signals 72 modeled by GPS car system module 26 drifted 
by as much as the equivalent of a i 1500 meter position error 
during signal loss from a particular satellite resulting from, 
for example obscuration in an urban setting, at. least one of the 
plurality of exclusive OR correlators 74 would provide the 
required prompt correlation to immediately lock onto satellite 
signals 72. 

Once the correlations are performed, the correlation results 
for each set of n Code Samples 80 are summed in summers 84 to 
produce a series of values each separately indicating the 
correlation of n Signal Samples 75 with each of the sets of n 
Code Samples 80. These correlation results are applied to 
threshold test 82, the output of which is applied to SatProcessor 
46 only when satellite signals 72 have been successfully 
received. The output of threshold test 82 specifies the number 
of delays which represent the prompt correlation for the 
reacquired satellite signal. It is important to note that in 
accordance with the present invention, the satellite tracking and 
reacquisition modes are not separated functions but rather 
interact seamlessly. That is, by providing a substantially 
expanded capture window, the correlations used for tracking are 
also automatically useful for immediate reacquisition as long the 
capture window is sufficiently wide to include any position error 
accumulated during signal obscuration or other loss. 

Because the speed of reacquisition is very important in 
order to maximize the opportunity to utilize the brief time 
during travel through intersection 22 when satellite 19 may 
temporarily be in view, it is advantageous to perform all such 
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correlations in parallel. Further, it is advantageous to 
continuously perform all such correlations in the capture window 
in order to minimize time when a satellite signal is not being 
tracked. In accordance with the presently preferred embodiment, 
exclusive OR correlators 74 are implemented in hardware rather 
than software to maximize the speed of correlation and minimize 
any error accumulation by minimizing the time for reacquisition . 
In operation, when car 10 follows track 62 along roadway 12, 
. during at least part of the time buildings 21 obscure the line of 
sight between car 10 and satellite 19. If satellite 19 had 
previously been acquired by GPS car system module 26, an 
approximate time value to synchronize with the satellite signals 
will be predicted. This value is maintained as accurately as 
possible within GPS car system module 26 while satellite 19 is 
obscured. In order to maintain the prediction for the required 
delay as accurately as possible, that is, to minimize the 
position error accumulated during signal loss, the above 
described techniques for maintaining or updating position 
accuracy by using cross-track hold, resetting position at a 
determined turn and/or the use of external sensors for dead 
reckoning provide a substantial benefit for use with the 
combined, expanded tracking and reacquisition windows described 
above . 

Present technology makes it convenient to provide H chip 
delays between correlators, but other delay values may be used. 
Similarly, it is convenient to expect that the prompt correlation 
can be maintained within plus or minus 5 chips of the timing of 
the satellite signals. Fig. 3 therefore portrays a series of 9 
or 10 early and 9 or 10 late correlators surrounding prompt 
correlator 74 to achieve the ± 5 chip capture window surrounding 
prompt correlator 74 in 20 half chip steps. A different number 
of correlators and other delays would also work with the present 
invention . 

Use of a plurality of fixed delays of one half chip width 
permit the immediate reacquisition of signals from a satellite to 
within an accuracy of one half chip width. In accordance with 
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satellite signals 72 as presently provided by the NAVSTAR 
satellites, one half chip width represents about 150 meters of 
maximum position error. It is possible to substantially reduce 
the maximum position error, and/or the speed of processing the 
data, by using fixed delays of a different amount of delay, e.g. 
fixed delays of one third, one quarter, one fifth or some other 
value of a chip width. 

Conventional approaches for different modes of operation, 
switch between wide and narrow delays at acquisition and/or 
reacquisition in order to provide a compromise between the width 
of the capture window and the number of correlations required for 
the desired range. In accordance with the present invention, a 
new technique is used which permits the convenient use of fixed, 
chip width delays to provide a finer gradation of correlation 
steps. In particular, as shown in Fig. 4, two sets of half width 
delays are used to provide the equivalent of a set of quarter 
width delays. The number of sets of fixed delays and the offset 
between them may be selected in accordance with the requirements 
of the application being addressed. 

Referring now to Fig. 4, a first plurality of sets of n Code 
Samples 80 are derived directly from code generator 7 6, delayed 
from each other by H chip width delays 78 and correlated with n 
Signal Samples 75 in exclusive OR (or NOR) correlators 74 as 
provided in Fig. 3. For convenience of explanation and drawing, 
the outputs from this first set of set of n Code Samples 80 are 
shown applied to summers 84 to indicate that the correlation 
products produced in exclusive OR correlators 74 from each such 
set of n Code Samples 80 are applied to threshold test 82 via 
summers 84. All such correlation products are applied, but for 
clarity only the correlation products having no delay, the 
predicted prompt or k/2th delay and the kth delay are depicted. 
The correlation products from this first plurality of sets of n 
Code Samples 80 are spaced apart by H chip width delay as noted 
above . 

In addition, in accordance with the present invention, 
additional sets of correlation products at different spacings are 
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available by use of one or more additional sets of ^ chip delays 
78 by, for example, tracking the same satellite in two or more 
channels offset in time from each other. It is important to note 
again that other delays and/or offsets may also conveniently be 
5 used and the delays need not all be the same., 

In particular, a second plurality of sets of n Code Samples 
are derived from code generator 76 and delayed from each other 
by H chip width delays 78. However, the delays in the second 
O sets of n Code Samples %A are offset from the delays in the first 

10 sets of n Code Samples 80 by a fixed amount, such as a % chip 
width delay, by insertion of H chip width delay 79 between code 
generator 76 and the first set of n code samples in sets of n 
d l5 Code Samples I4f. t This results in each of the samples in sets of 
£ ;Ej n Code Samples %X falling halfway between two of the sets of n 

ijijlS Code Samples 80. As shown in Fig. 4 only k-1 sets of n Code 
(I N[i Samples are required with k sets of n Code Samples 80. 

; p Each of the sets of n Code Samples 84 are correlated with n 

5j3 Signal Samples 75 in exclusive OR correlators 74 as provided in 
Fig. 3 to produce correlation products which are then summed by 
ny20 additional summers 84. As noted above, the dashed lines between 
^ each of the sets of code samples and summers 84 are used to 
\m indicate that the correlation product between that set of code 
13 samples and n Signal Samples 75 is applied to a particular one of 
summers 84. As can then easily be understood from Fig. 4, 
25 correlation products separated from each other by H chip width 

delays, from the 0th delay to kth delay, are produced using sets 
of H chip width delays and a single H chip delay (which may 
represent the offset delay between two channels) and after 
individual summation are applied to threshold test 82 to 
30 determine which delay represents the currently prompt delay of 

satellite signals 72 from a satellite being reacquired by GPS car 
system module 26. 

The second set of H chip delays may easily be implemented by 
having a second channel track the same satellite, offset, however 
35 by H chip width delay 79. 

In this way, the range of delay within which a satellite 
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signal lock may be acquired, maintained and/or reacquired may be 
reduced from ± H chip width, to about ± H chip width, which 
permits faster pull in to lock, i.e. when the tracking has been 
optimized and range error reduced to minimum. 

It is important to note the seamless integration of tracking 
and reacquisition provided by the present invention in that the 
same correlations are used for tracking and reacquisition and the 
related speed of capture and lock and simplicity provided 
thereby. The ability to rapidly reacquire within a capture 
window so that one of the correlations may immediately be used as 
a prompt correlation, speeds up all data acquisitions thereafter. 

It is also convenient to utilize a first plurality of sets of n 
Code Samples 80 for tracking and, when satellite signals 72 are 
lost, provide additional accuracy in reacquisition by using a 
second plurality of sets of n Code Samples such as sets of n Code 
Samples 84. In particular, the same plurality of sets of n Code- 
Samples 84 may be used for reacquisition of signals 72 for 
different satellites at different times in order to reduce the 
total number of components and steps required to produce all the 
necessary correlations and summations. 

In operation, GPS car system module 26 continuously attempts 
to track and reacquire the signals from satellite 19 in SatTRAK 
channel 38 while satellite 19 is obscured from view. As car 10 
passes through intersection 22, the line of sight to satellite 19 
is momentarily not obscured by buildings 21. Whenever any of the 
correlations performed in SatTRAK channel 38 indicate . that the 
satellite signals are being received with sufficient strength so 
that the correlation products from some of the correlators are 
above threshold, reacquisition is immediately accomplished. 
Reacquisition occurs when the correlator output indicating the 
largest magnitude is selected as the new prompt correlator. 
Conventional techniques for improving the quality of the data are 
then employed. 

The data from satellite 19 is used to immediately, after a 
settling time for lock, update the GPS data and correct the 
currently known position information derived therefrom. Even if 



satellite 19 is then again immediately obscured, the update 
information derived during travel through the intersection by 
fast reacquisition provides a substantial improvement in accuracy 
of the GPS determined position. This will permit GPS car system 
module 2 6 to continue accurate navigation even through otherwise 
very difficult areas, such as city streets. 

Although the use of single satellite navigation data by 
cross-track hold and then updating a satellite data by detecting 
turns and/or immediately reacquiring satellite signals in 
intersections have all been described separately, they are also 
very useful in combination. Terrestrial navigation systems, 
using GPS receivers in a stand alone mode, aided by map displays 
and data bases and/or aided by external sensors such as inertial 
navigation systems may benefit from the use of combinations of 
one or more such modes. In a preferred embodiment of the present 
invention, all three techniques are combined to maximize the 
ability of the car navigation system to provide accurate and 
useful navigation data while traversing a difficult environment 
such as city streets. 

Referring now to Fig. 5, a preferred embodiment of the 
present invention is described in which major portions of SatTRAK 
channels 38, 40, 42 and 44 and SatProcessor 46 of the present 
invention are implemented in an Application Specific Integrated 
Circuit or ASIC 102. Many of the functions of a conventional 
satellite processor may still, however, be performed in software. 
The particular implementation depicted provides a 12 channel GPS 
acquisition and tracking system with fast reacquisition 
capabilities as described above while substantially reducing the 
number of gates required on the ASIC to implement this system. 

The signals received by GPS antenna 28 are digitized and 
downconverted to form a digital composite of signals received 
from all satellites in view to produce sample data 100 which is 
at a frequency of 37.33 f 0 where f 0 is the chip rate of the C/A 
code applied to each GPS satellite. For convenience, the 
frequencies described below will be designated in terms of 
multiples of f 0 . Each of 12 Space Vehicles (SVs) or satellites 
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are tracked in ASIC 102 under the control of Central Processing 
Unit, or CPU, 101 which provides control signals and data to ASIC 
102. In particular, CPU 101 provides data regarding the 
predicted Doppler shifts and C/A code applied to each SV to 
5 Random Access Memory, or RAM, Rl 103 associated with ASIC 102 
which provides the data to RAM R2 105 at designated times. RAM 
R2 105 provides data to and receives data from ASIC 102, 
permitting CPU 101 data updating and ASIC 102 processing of old 
data to operate simultaneously. RAM R2 105 is used as a buffer 

10 by ASIC 102 primarily to store intermediate values of signals 
during processing. Other conventional portions of a micro- 
computer including a CPU are not shown but conveniently may 

5 j3 include devices operating software implementing the single 

satellite, cross-track hold and other techniques described above 

& as well as other functions of SatProcessor 46. 

ny Sample data 100 is applied to C/A code acquisition, tracking 

and reacquisition block CACAPT 104 in ASIC 102 where it is split 
into in-phase and quadrature-phase, or I and Q, signals at 
j]r baseband by I/Q splitter 106. After processing by CACAPT 104, 
gp the I,Q signals are rotated for Doppler shift in 12 channel 
i*= Doppler Block 108 which separately compensates for the expected 
!S Doppler frequency shifts of each of the 12 SVs which can be 
tracked. 

The Doppler rotated I,Q signals for each SV are then applied 
25 to Correlator Block 110 where each signal sample, which is from 

one of the 12 SVs, is correlated in a multiplexed fashion with 20 
delayed versions of the C/A code, produced by 12 channel Coder 
Block 112, for that SV. During each segment of time, as 
described below with regard to Fig. 11 in greater detail, 
30 Correlator Block 110 performs 240 C/A code correlations in 
accumulator 175 to enhance the speed of acquisition and 
reacquisition. The output of Correlator Block 110 is applied to 
IQACCUM Block 114, and the output of IQACCUM Block 114 is applied 
to IQSQACCUM 116, in Accumulator Block 115. IQACCUM Block 114 is 
35 conveniently configured from another block of RAM associated with 
ASIC 102, identified herein as RAM 3. Similarly, IQSQACCUM 116 
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is conveniently configured from another block of RAM associated 
with ASIC 102, identified herein as RAM 4. 

Accumulator Block 115 operates in different fashions during 
acquisition, tracking and reacquisition modes under the direction 
of CPU 101. During acquisition mode, Coder Block 112 is caused 
to sequence through as many sets of 240 different codes delays as 
necessary to acquire the satellite signals from a particular 
space vehicle. That is, as many sets of 240 different delays are 
correlated in Correlator Block 110 to provide IQSQACCUM 116 with 
an appropriate correlation output whose power indicates that 
correlation has been achieved with that satellite. The process 
is then repeated for each satellite to be acquired. For 
convenience, all delays may be tested. 

During reacquisition, a single set of 20 delays are 
correlated in Correlator Block 110 to determine if one such delay 
provides a peak value above a predetermined threshold to indicate 
that a correlation has been achieved and the satellite thereby 
reacquired. The reacquisition mode operates transparently within 
the tracking mode in that a set of 20 delays are correlated in 
Correlator Block 110. If tracking is maintained, the peak signal 
may migrate from a particular delay to the next adjacent delay 
but will be maintained within the current set of 20 delays being 
correlated. It is convenient to consider the delay producing the 
signal with the greatest magnitude as the prompt delay, producing 
the prompt correlation product. The signals produced by one more 
and one less delay then become the early and late correlation 
products which may be processed in a conventional manner to 
maintain lock with each satellite. 

If the- signal from the satellite is temporarily obscured or 
lost for any other reason, the then current set of 20 delays is 
correlated and searched for a peak of sufficient magnitude to 
indicate reacquisition. The Doppler and code values are 
continuously updated based upon the last available position 
information including velocity, and the correlations are 
performed, until the satellite signal is reacquired or sufficient 
time has elapsed so that the satellite signal is considered lost. 
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The operation and configuration of ASIC 102 will now be 
described in greater detail with regard to the in-phase or I 
signal path. The quadrature phase or Q signal path is identical 
and need not be separately described. 

Within CACAPT 104, sample data 100 is applied at 37.33f 0 to 
I/Q splitter 106 to produce a 2 bit signal at 18.67f 0 which is 
further reduced to 2f 0 by Digital Filter 118 which operates by 
adding sets of 10, 9 and 9 samples which are summed, quantized, 
and then stored serially in 11 sample deep buffer 120. When 11 
sample deep buffer 120 is filled, the data is transferred in a 
parallel fashion to an identical buffer, called parallel block 
122, for Doppler rotation. Data is therefore transferred out of 
11 sample deep buffer 120 when 11 samples are received, that is, 
at a chip rate of 1/llth of 2f 0 or approximately 0.18f o . 11 
sample deep buffer 120 operates as a serial to parallel converter 
while parallel block 122 operates as a parallel to serial 
converter. This results in 186 parallel transfers per msec. 

Data is shifted out of parallel block 122 at 24f 0 to 12 
channel Doppler Block 108 so that the Least Significant Bit- or 
LSB of the serial converter, parallel block 122, is the output of 
CACAPT 104 in the form of CapIOut and CapQOut which are applied 
as CACAPT Data output 123 to 12 channel Doppler Block 108. The 
increase in chip rate from 2f 0 to 24f 0 provides an operating 
speed magnification of 12 as will be described below in greater 
detail . 

Referring now also to Fig. 6, 12 channel Doppler Block 108 
is now described in greater detail. Doppler Block 108 receives 
satellite specific CACAPT Data output 123 including CapIOut and 
CapQOut from CACAPT 104 for storage in Doppler Register 124. 
Satellite or source specific predicted Doppler phase, after 
processing by Carrier Numerical Control Oscillator or NCO 125 and 
sine/cosine look-up table 134, is also applied to Doppler 
Register 124 where it is added to CapIOut and CapQOut for the 
same SV (or other source) to form dopIOut and dopQOut . Within 
Doppler Block 108, Carrier_NCO 125 operates at an effective rate 
of 2f 0 for each satellite channel because the data sample rate is 
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For each SV, CPU 101 stores the satellite specific predicted 
carrier phase dopPhaseParam, and predicted carrier frequency 
dopFreqParam, in RAM R2 105. Sat_Mem 186 then transfers the 
dopPhaseParam and dopFreqParam as shown in Fig. 9 to Carrier 
Phase Register 126 and Carrier Phase Output Buffer 128, 
respectively, at each 1 msec boundary. Adder 130 adds carrier 
phase to carrier frequency, derived from dopPhaseParam and 
dopFreqParam, to produce the current carrier phase value in 
Carrier Phase Register 126 shown as Carrier_NCO. 

The four Most Significant Bits or MSBs of Carrier_NCO in 
Carrier Phase Register 126 are applied to sine/cosine look-up 
table 134 which includes 2 4-bit registers for storing its 
output. The output of sine/cosine look-up table 134 is applied 
to Doppler Multiplier 132 in Doppler Register 124 for Doppler 
rotation of CACAPT Data output 123 (CapIOut and CapQOut) to 
produce rotated SV output signals dopIOut and dopQOut. Doppler 
Register 124 uses Doppler Multiplier 132, as well as four 4-bit 
registers, two adders, another pair of 5-bit registers and a 
quantizer to form dopIOut and dopQOut . Referring for a moment to 
Fig. 8, dopIOut and dopQOut are applied to parallel converter 166 
and rotated SV output signal 127 is the output of serial to 
parallel converter 166 which is applied directly to 11 bit 
Holding Register 140. 

During each segment of time, the beginning value for the 
Doppler phase of each SV is -stored in RAM R2 105, retrieved 
therefrom by Doppler Block 108 for the rotation of the SV during 
that segment. At the end of each segment, the end value of 
Doppler phase is stored in RAM R2 105 for use as the beginning 
value for the next segment. Under the control of gpsCtl 182, 
Doppler phase value dopP_Next in Carrier Phase Output Buffer 128, 
saved at the end of each rotation for a particular SV by dopSave, 
is applied to Sat_Mem 186 for storage in RAM R2 105 for that SV, 
to be retrieved by Doppler Block 108 again during the next 
Doppler rotation of that SV in the following segment. The 
operation of Multiplexer Block 129 may be best understood from 
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the description of the triple multiplexing of ASIC 102 associated 
with Fig.s 10 and 11. 

Referring now also to Fig. 7, 12 channel Coder Block 112 
includes Coder_NCO 136 and Code Generator 138. Coder_NCO 136, 
which is similar to Carrier_NCO 125 shown in Fig. 6, creates 
Gen_Enable whenever Phase Accumulator 148 overflows. Gen Enable 
is the MSB of the output of Phase Accumulator 148 and is applied 
to Code Generator 138. 

In particular, under the control of gpsCtl 182, Sat_Mem 186 
applies the satellite specific 24 bit code frequency parameter, 
coderFreqParam, and the 24 bit satellite specific code phase 
parameter, codePhaseParam, at each 1 msec edge to Coder NCO 136 
from RAM R2 105. CoderFreqParam is added to codePhaseParam 
effectively at 4f 0 per channel in Phase Adder 150 even though 
codePhaseParam operates at 48f 0 during tracking and 
reacquisition. A pulse can be generated for Gen_Enable between 0 
Hz and 4f 0 Hz. In order to generate Gen_Enable at 2f 0 , the value 
of half the bits (23:0) of Phase Accumulator 148 must be loaded 
in as coderFreqParam. 

The LSB of codePhaseParam represents l/256th of a C/A code 
chip. CodePhaseParam initializes the contents of Phase 
Accumulator 148. Gen_Enable is generated whenever Phase 
Accumulator 148 overflows. Phase Accumulator. 14 8 is a 25 bit 
register initialized by the value of codePhaseParam when 
corHoldRegLoad 152 from CPU 101 is active at each 1 msec edge 
when new data is written from CPU 101, The 24 LSBs of 25-bit 
Phase. Accumulator 148 are then added to coderFreqParam in Phase 
Adder 150 and returned to Phase Accumulator 148. Phase Buffer 
Register 154 stores and buffers the contents of Phase Accumulator 
148, to produce CoderPNext which is updated whenever codCodeSave 
158 from gpsCtl 182 is active. CoderPNext is applied to Sat_Mem 
186 for storage in RAM R2 105. The operation of multiplexer 142 
may be best understood from the description below of the triple 
multiplexing of ASIC 102 provided with Fig.s 10 and 11. 

Gen_Enable is applied to Code Generator 138 to cause a new 
• code to be generated. C/A Codes parameters Gl and G2 are 
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parallel loaded from RAM R2 105 by Sat_Mem 186 as glParln and 
g2ParIn into Code Generator 138 to produce glGenOut and g2GenOut 
which are returned to RAM R2 105 by Sat__Mem 186. The bit-0 of 
both Gl and G2 generators in Code Generator 138 are internally 
5 XOR=d and generate genSerOut 160 which is serially applied to 11 
bit Code Shift Register 170 in Correlator Block 110, as shown in 
Fig. 5. Code Generator 138 generates the following C/A codes: 
Gl = 1 + X3 + X10 

G2 - 1 + X2 + X3 + X6 + X8 + X9 + X10. 
10 The output of Code Shift Register 170 is applied to correlators 

74, 11 bits at a time at 48f 0 so that at least 20 code delays, 
^ separated by one half chip width, are correlated against each 
:jfj Doppler rotated sample from each SV. The increase in chip rate 
^ from 2f 0 to 48f 0 provides a magnification factor of 24 as will be 
ill described below in greater detail. 

nil Values of Gl and G2 are be stored in RAM R2 105 during each 

^ segment after correlation with the Doppler rotated sample in 

correlators 74 for that SV so that they may then be retrieved by 
^ Coder Block 112 during the next time segment for correlation of 
^tj the next 11 bit sample from the same SV. 

Referring now also to Fig. 8, Correlator Block 110 is shown 
!S in greater detail. . DopIOut and dopQOut in the rotated SV output 
from Doppler Block 108 are applied to serial to parallel 
converter 166 which is then parallel loaded to Holding Register 
25 140. GenSerOut 160 from Coder Block 112 is applied to Code Shift 
Register 170 in Correlator Block 110. These data sets represent 
the Doppler shifted data received from the SV, as well as the 
locally generated code for that SV, and are applied to Exclusive 
NOR gate correlator 74 for correlation under control of gpsCtl 
30 182. 

The output of correlator 74 is applied to Adder 174 and 
combined in Bit Combiner 176 to corlOut 178 and corQOut 180 which 
are applied to IQACCUM Block 114 and IQSQACCUM 116 shown in Fig. 
5. Adder 174 and Bit Combiner 176 operate as a partial 
35 accumulator as indicated by accumulator 175 in Fig. 5. 

Referring now also to Fig. 9, on overview of the operation 
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of ASIC 102 is shown. A dedicated set of on-chip logic controls 
the operation of ASIC 102 and is identified herein as gpsCtl 182 
In particular, under the control of gpsCtl 182, sample data 100 
from the GPS satellites is applied to CACAPT 104 where it is 
separated and decimated into I and Q data streams to form CACAPT 
Data output 123. SV data 123 is rotated for the predicted 
Doppler shift of each SV to produce rotated SV output signals 
dopIOut and dopQOut which are correlated with genSerOut 160 from 
Coder Block 112 in correlators 74. CorlOut 178 and corQOut 180 
from correlators 74 are accumulated in IQACCUM Block 114 and 
IQSQACCUM 116 to produce output 184 to CPU 101. 

As will be further described below in greater detail, a 
portion of memory is used for SatJMem 186 which stores and 
provides the Doppler shift and code information required during 
multiplexing. 

In operation, every millisecond is divided into 186 
segments, each of which includes 264 clocks. Within these 264 
clocks, 12 channels are processed with each channel taking 22 
clocks to compute 22 different correlations or delays. Only 20 
of these 22 correlations are stored and used for subsequent 
processing. For each channel, gpsCtl 182 controls the loading of 
Carrier_NCO 125 in Doppler Block 108 using dopLoad and dopSave. 
Similarly, gpsCtl 182 controls the loading of Coder_NCO 136 in 
Coder Block 112 via corHoldRegLoad and corCodeSave. The flow of 
data through Correlator Block 110 is controlled with 
serialShiftClk, and also corHoldRegLoad and codCodeSave. Control 
signals are applied to IQACCUM Block 114 and IQSQACCUM 116 for 
each channel and include startSegment, startChan, resetAcc, peak, 
iqsq, wrchan, ShiftSellqSq and acq_mode. Within each segment, 
gpsCtl 182 provides the periodic signals eng_capShif tClk, 
capLoad, syncpulse, serialShiftClk to CACAPT 104 to repackage 
incoming satellite data samples into groups of 11 half chip 
samples. 

All accesses initiated by gpsCtl 182 are processed by 
SatJYfem 186 to generate read/write control and address signals 
for RAM Rl 103 and RAM R2 105. GpsCtl 182 controls the flow of 
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data through all data paths together with Sat_Mem 186 and manages 
the access of channel parameters stored in RAM Rl 103 and RAM R2 
105. RAM Rl 103 is written to by the user to define the channel 
parameters that will be loaded to RAM R2 105 at the end of the 
corresponding integration or accumulation time. RAM R2 105 is 
used by the data path as a scratchpad to store the intermediate 
values of the various channel parameters during processing. 

Data read out of RAM R2 105 is sent to the various parameter 
registers in Doppler Block 108, Coder Block 112, Correlator Block 
110 and gpsCtl 182 under the control of Sat_Mem 186. Data from 
these blocks and RAMI 190 are multiplexed at the input to the 
write port of RAM R2 105. RAM Rl 103 is a 16x108 asynchronous 
dual port ram used for the parameters for all 12 channels while 
RAM2 192 is another 16x108 asynchronous dual port ram used for 
storing intermediate values of the satellite parameters during 
processing, while switching from one channel to the next. 

Referring now to Fig. 10, the system of the present 
invention includes a multiplexed data path in order to reduce the 
size and complexity of ASIC 102 on which the majority of the 
parts of the system can be provided. Conventional receiver 
designs have multiplexed a single set of correlators for use for 
each of the separate channels in which. an SV is tracked in order 
to reduce the number of correlators required. The use of the 
system of the present invention reduces the million or more gates 
that would be required for a conventional configuration down to a 
manageable number, on the order of about less than 100,000. 

In accordance with the present invention, in addition to 
multiplexing the satellite channels in a manner in which no data 
is lost/ the code delay correlations are also multiplexed. That 
is, conventional receivers use two or three correlators to 
provide early, late and/or prompt correlations for each SV. The 
present invention multiplexes a plurality of code delays in order 
to provide far more code delay correlations than have been 
available in conventional systems without substantially 
multiplying the hardware, or chip area on ASIC 102 required by 
the number of gates used. 



33 



10 



The multiplexing of code delays permits the wide capture 
window described above with regard to Fig.s 3 and 4 that permits 
rapid SV reacquisition . In particular, 20 delays such as *s chip 
delays are provided and constantly monitored for each SV so that 
GPS data can be acquired even during brief glimpses of the SV, 
for example, when car 10 is in intersection 22 as shown in Fig. 
1. The SV can be reacquired and useful data obtained because the 
modeling of the vehicle's position on roadway 12 is sufficiently 
accurate to keep the predicted code and Doppler values for a 
previous acquired and currently obscured SV within a window of ± 
10 half chip code delays. In this way, data obtained during 
>h reacquisition can be used directly as GPS data. That is, the 
,£ reacquisition mode is transparent to the tracking mode. The GPS 
data is acquired whenever available without substantial lost time 
for reacquisition. 
Hy! Further, the operation of satellite tracking is itself 

;| multiplexed for each set of data for all 12 channels in order to 
further substantially reduce the ASIC gate count. That is, only 
!* a small portion of the bits in the C/A code is processed at one 
2j time for all 12 SVs . In order to digitally process the signals 
i* received, the digital representations of these signals must be 
;|" processed in registers and buffers capable of storing the digital 
data. The C/A code contains 1023 bits in each repetition which 
lasts 1 msec. If all 1023 bits were to be processed at once, 
25 registers 1023 bits wide would be required. Such registers would 
be expensive in cost and gate count and quite cumbersome. In 
accordance with the third level of multiplexing used in the 
triply multiplexed receiver configuration of the present 
invention, a smaller register is multiplexed to handle different 
30 portions of the 1023 bits of the C/A code. This means the 

smaller register is used many times during each 1 msec repetition 
of the C/A code to process enough smaller samples of the data 
received so that within each msec all 1023 bits can be processed. 
In the preferred embodiment described above particularly in 
35 Fig.s 3 to 9, a configuration using 11 bits registers was used so 
that each register is used 186 times per msec to process all 1023 
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bits of a C/A code repetition. Each l/186th of a msec is called 
a segment. The tracking of each SV is therefore multiplexed 186 
times by processing the 11 bits in each register during each 
segment. In addition, in the preferred embodiment, 12 channels 
5 are used to track a maximum of 12 SVs. This requires that each 
11 bit segment is multiplexed 12 times during that segment to 
apply a Doppler rotation for each SV. 

Further, each channel is further multiplexed by a factor of 
22 to provide a substantial plurality of different code delays. 

10 This requires that the Doppler rotated sample for each SV is 
correlated 22 times with different C/A Code delays before the 
Doppler rotated sample for the next channel is produced. In this 

^ manner, 22 different code phases may be tested for each of 12 SV 
during each of 186 segments to provide real time data with only 

0|S 11 bit wide registers by processing each register 186 times per 

ill msec. 

: f It is important to note that the processing of the present 

'~ invention occurs during a particular segment, i.e. a l/186th of a 
1^ repetition of the C/A code, during the length of time required 
2jb for the segment to be collected. In this optimized manner, no ' 
1^ data is lost during tracking or reacquisition or switching 
ig between these states because the data being processed in any 
particular segment is at most 11 half chip delays old. 

Referring now to Fig.s 10 and 11, the output of Digital 
25 Filter 118 shown in Fig. 5 is sample data stream 119 at 2f 0 . The 
chip rate of the C/A modulation of the signals 100 from the SVs 
is at f 0 . In order to avoid loss of any data, the SV signals 
must be sampled at least at their Nyquist rate, that is, at twice 
the chip rate of the modulation of interest which is 2f 0 . 
30 Although sample data stream 119 can be operated at a higher chip 
rate than the Nyquist rate, which is twice the chip rate, there 
is no advantage in doing so. 

Sample data stream 119 is therefore a series of samples of 
the digitized and filtered SV data at twice the chip rate of the 
35 C/A code, that is, each sample in sample data stream 119 has a 
width equal to one half of a C/A code chip. The number of bits 
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in each msec or cycle of code in sample data stream 119 is twice 
the number of bits in the modulation, i.e. 2046 bits each 
representing one half of a C/A code chip. In accordance with the 
multiplexing scheme of the preferred embodiment being disclosed, 
the data is processed in 11 bit segments, and sample data stream 
119 is therefore applied serially to 11 bit (10:0) register value 
buffer 120. The time required to serially store 11 bits out of a 
total of 2046 bits in the 2f 0 data stream is 1 + (2046 + 11 = 
186) or l/186th of a msec. 

During the time the first set of 11 sample bits is being 
stored in 11 sample deep buffer 120, no bits are available for 
processing. After the first 11 sample bits are serially received 
and serially stored, the 11 sample bits are transferred in 
parallel to parallel block 122. This parallel operation 
therefore occurs every l/186th of a msec or at a rate of 
approximately 0.18f 0 . Each l/186th of a msec is called a time 
segment or segment and is the unit of processing for most of the 
operations. The 1023 chip C/A code of each of the satellites in 
the composite signal received is processed in 11 half chip bits. 

Dividing the msec repetition rate of the C/A code into 186 time 
segments multiplexes each of the 11 bit registers by a 
multiplexing factor of 186. 

CACAPT Data output 123 from parallel block 122 is processed 
in Doppler Block 108 at a much faster chip rate, for example at 
24f 0 . That is, the 11 bits of sample data in each segment of 
time is multiplexed by a factor of 12 to permit 12 different 
operations to be performed to that set of 11 bits of data. In 
particular, in Doppler Block 108, CapIOut and CapQOut of CACAPT 
Data output 123 are multiplied in Doppler Register 124 by twelve 
different Doppler shifts so that within each segment twelve 
different Doppler rotations are performed. 

Each different Doppler shift represents the predicted 
Doppler rotation required for each of the maximum of 12 different 
SVs that can be tracked. The increase in processing chip rate 
from 2f 0 to 24f 0 multiplexes the processing for each of 12 
channels of data. It is important to note that the multiplexing 



to permit one channel to operate as 12 multiplexed or virtual 
channels each representing a different SV is applied only after 
the input signals are multiplexed, that is, broken into 186 time 
segments each including 11 half chip width bits. In this way, 
the multiplexing for 12 channels or satellites is easily 
accomplished with relatively inexpensive 11 bit registers without 
loss of time or data. The selection of the number of sampling to 
be an integer division of the number of code bits per period is 
important to achieve these goals. Multiplexer Block 129 in 
Carrier_NCO 125 controls the timing of this multiplexing under 
the direction of gpsCtl 182. 

The output of Doppler Block 108, signals dopIOut and 
dopQOut, are applied to serial to parallel converter 166 within 
Correlator Block 110. Each rotated SV output signal 127 
represents the rotated signal from a single SV and 12 such 
rotated SV output signals 127 are produced in each segment of 
time • 

. Rotated SV output signal 127 is. loaded in parallel fashion 
into Holding Register 140 in Correlator Block 110. The input to 
Exclusive NOR gate correlator 74 is therefore an 11 bit wide 
signal which is retained for l/12th of a time segment as one 
input to Exclusive NOR gate correlator 74. 

Correlator 74 is a series of 11 separate one bit correlators 
which all operate in parallel. One input is rotated SV output 
signal 127 while the other 11 bit input is provided by 11 one bit 
genSerOut 160 output bits from Coder Block 112. During the 1/12 
of a time segment provided for operation on the rotated SV output 
signal 127 for a particular satellite, the code for that SV is 
produced serially by Code Generator 138 and applied to Code Shift 
Register 170. 

At the beginning of the correlation for a particular 
channel, 11 bits of the code for that SV have been shifted into 
Code Shift Register 170 and are available therein for 
correlation. Every l/22nd of a channel (that is, a 1/12 of a 
segment) each of the 11 bits in Code Shift Register 170 are 
correlated in one of 11 one bit exclusive Nor gates in Exclusive 



NOR gate correlator 74. This produces 11 correlator output bits, 
the sum of which indicates the magnitude of the correlation 
between the rotated SV output signal 127 and that code phase. 
These 11 correlation sums produced in parallel are summed in 
parallel and stored in the first of 22 summers related to that SV 
in Accumulator Block 115. 

During the next or second l/22nd of a channel, Code 
Generator 138 produces the next bit for the C/A code for that SV. 

This next bit is applied serially to Code Shift Register 170. 
At this time, 10 bits from the first correlation remain in Code 
Shift Register 170 and together with the newest bit form another 
11 bit sample of the expected code for that SV, delayed from the 
previous 11 bit sample by the time required to generate 1 bit, 
that is, one half chip width at the rate code is produced, 48f 0 . 

The second sample is therefore a one half chip delayed version 
of the code, delayed one half chip width from the previous 11 bit 
samples. It is important to note that the two 11 bit code 
samples just described differ only in that a new bit was shifted 
in at one end of the register to shift out the MSB at the other 
end of the register. 

The 11 bit correlation product of the same rotated SV output 
signal 127 and the second 11 bit sample of code is then stored in 
the second of the 22 summers related to that SV in Accumulator 
Block 115. Thereafter, the remaining 20 serial shifts of the 
genSerOut 160 from Code Generator 138 are correlated against the 
same rotated SV output signal 127 to produce 20 more sums of 11 
bit correlations for storage in Accumulator Block 115 for that' 
SV. The result is that 22 values are then available within 
Accumulator Block 115 for processing, each value is a measure of 
the correlation of the signals from one SV with 22 different code 
phases or delays, each separated by one half chip width. 

During the next 1/12 of a time segment, that is, during the 
processing of the second multiplexed channel, the rotated SV 
output signal 127 for the next SV, is applied to Holding Register 
140 for correlation with 22 different one half chip delays of the 
code generated for that satellite. At the end of a segment, 
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Accumulator Block 115 includes a matrix of 12 by 20 different 
sums. In one implementation of the present invention, it has 
been found to be convenient to save only 20 out of the 22 
possible code delay correlation results. The 12 rows of 20 sums 
represent the measure of correlation for each of the 12 SVs at 20 
code phases or delays. 

In summary, the data path for the present invention is 
triply multiplexed in that 

(a) each msec, which represents 1023 bits of C/A code, 
is sliced into 186 to form the 186 segments in a msec of sample 
so that only 11 half chip wide sample bits are processed at one 
time; 

(b) each segment is then multiplexed by 12 so that each 
such 11 bit sample is rotated for twelve different sources; 

(c) the rotated 11 bit. sample for each source is 
correlated against 20 sets of different code delays for that 
source to multiplex within each channel by 20; and 

(d) the sum of the correlation products for each delay 
in each channel are then summed to produce the. accumulated 
correlation output. 

Although 22 different delays are available, it is convenient 
to use 20 such delays, or code phase theories for testing the 
rotated satellite signal. The correlation product having the 
greatest magnitude for each channel after accumulation, that is, 
the largest of the 20 sums of 11 bits stored in Accumulator Block 
115 for each channel may then be detected by its magnitude, for 
example by a peak detector, to determine which delay theory is 
the most accurate. The peak sum represents the on-time or prompt 
correlation for that SV. 

Turning now specifically to Fig. 11, the triple multiplexing 
scheme of the present invention may easily be understood by 
looking at the slices of time resulting from each of the 
multiplexing operations. Within each msec, the C/A code for each 
particular satellite has 1023 bits. In order to preserve all 
necessary information, the satellite signals are sampled, in a 
digital composite of signals from all satellites, at the Nyquist 



rate at 2f 0 to produce 2046 half chip wide sample bits. 

Each sequential set of eleven sample bits are processed 
together as a segment of time, the length of which is equal to 1 
/ (2046 + 11) of a msec, i.e. one l/186th of a msec. After 
processing of the 186th segment in a msec all necessary data has 
been extracted and the 11 bit sample for the next segment is 
available. Although the partial sums accumulated over each msec 
in Accumulator Block 115 may only be evaluated at the end of a 
msec, no data is lost and the results are only 1 segment late. 
That is, since it takes 1 segment to fill 11 sample deep buffer 
120 and transfer the 11 bit sample to parallel block 122, the 
data from the first 11 bit sample is being processed while the 
data for the second 11 bit sample is being collected. Even if 
the system operated for a year, the sampled being processed to 
provide position information is still only one time segment old. 

The 11 bits of each segment are multiplexed for each SV by 
being time division multiplexed during Doppler rotation. That 
is, the 11 bit sample of segment 1 is used to provide 12 
different Doppler shifted outputs so that a single 11 bit segment 
sample is used 12 times to produce 12 different satellite 
specific Doppler rotated versions, assuming all 12 satellites are 
in view or being modeled. The operations for one channel then 
require one twelfth of a segment. It is critical to note that 
each segment only produces a partial result and that the 12 
partial results during each segment must be summed at the end of 
each msec to provide valid output data. 

Each of the operations on one particular channel in a 
segment are time division multiplexed by a factor of 22 so that 
22 different code delays for that partial sum for that satellite 
can be tested. The peak sum of these 22 correlations can however 
be detected by magnitude immediately if necessary to select the 
most likely delay for that channel. In the present embodiment, 
the information for that channel is only valid once per msec when 
summed or accumulated so that there may not be a substantial 
advantage in peak detected with a particular segment. In some 
GPS applications and in other spread spectrum applications, such 
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as' wireless communications, it may be desirable if strong signals 
are present to accumulate and transfer the sum of the 
accumulations for each source from R3 to R4 more often than once 
per code repetition rate. The time required to evaluate a 
5 particular code phase delay or theory is only l/22nd of the time 
required per channel per segment or l/22nd of 1/12 of l/186th of 
a msec. This speed of operation is more easily achieved because 
the 11 one bit correlations required are produced in parallel. 
Similarly, the speed of generation of the different code delays 
10 for a particular SV is more easily accomplished in accordance 

with the present invention because each 11 bit code delay sample 
is automatically produced when each single new bit, i.e. each new 
:j5 genSerOut 160, is shifted into Code Shift Register 170. 
* s% = The selection of the magnitudes or multiplexing factors used 

X|* in each level of multiplexing is not arbitrary. The larger the 
jiy number of segments, the smaller the required size or depth of the 
registers need for each sample. By using a code repetition 
multiplexing factor of 186, that is, by dividing the 2046 bits of 
a 2f 0 by 186, only 11 sample bit, need to be evaluated at a time. 
2;(J The number of required channels is bounded pragmatically by 

h= the fact that at least 4 SVs must be in view at the same time to 
!;£ determine position accurately in three dimensions. Time is the 
fourth unknown which must be determined along with each of the 
three dimensions although provisions for estimating, modeling 
25 and/or updating the position information as described above so 

that position information may be accurately provided even during 
periods when less than 4 satellites are concurrently in view. 

The constellation of 24 NAVSTAR . satellites in use are 
arranged to cover the earth so that a maximum of 12 such 
30 satellites may be in view at any one location at any particular 
time. The maximum number of pragmatically useful channels is, 
for this reason, no less than about 12 channels. The selected 
channel multiplexing factor used in the channel level of 
multiplexing in the embodiment shown herein is therefore a factor 
35 of 12. 

The number of different code delays is bounded at the low 
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end by an absolute minimum of 1 so that if the exact delay can 
somehow be maintained, the only necessary correlation would be 
the on-time or prompt correlation. Conventional GPS receiver 
systems use at least 2 or 3 different code delays so that 
5 conventional tracking techniques, for example those which use 
early, prompt and late correlations to center the prompt 
correlation within ± 1 delay, may be employed. 

In accordance with the present invention, a substantially 
greater number of different code delays, or delay theories, are 
10 tested so that fast reacquisition may be accomplished as 

described above with regard to Figs. 3 and 4. Although for the 
. 3 H particular preferred embodiment described herein, it was 
ifl determined that a total of 20 different delays, each separated in 
^ time by one half the width of a C/A code chip, i.e. 1/2 of 1/204 6 
lifj of one msec, a code delay multiplexing factor of 22 was selected 
ny because the relationship between each of the 3 multiplexing 
; |h factors is also important. 

The product of the three multiplexing factors, code 
:]n repetition multiplexing factor, channel multiplexing factor and 
2;(i code delay multiplexing factor should optimally be an even 
j J integer multiple of the number of bits in each repetition of the 
spread spectrum modulation. An even integer multiple is required 
because samples must be taken at twice the chip rate, i.e. at the 
Nyquist rate, in order to avoid data loss from sampling at a 
25 slower rate. Although multiplexing factors can be used 

successfully even if the product is not exactly equal to an even 
integer multiple, data- loss or unnecessary complexity and costs 
may result. 

In the particular embodiment shown, the spread spectrum code 
30 of interest is the C/A code, each repetition of which includes 
1023 bits. In accordance with the triple multiplexing product 
rule discussed above, the product of the three multiplexing 
factors must equal an even integer multiple of 1023, such as 
2046. In the described embodiment, the code repetition 
35 multiplexing factor is 186, the channel multiplexing factor is 12 
and the code delay multiplexing factor is 22. The product of 186 
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multiplied by 12 and then by 22 is 49104 which, when divided by 
1023, equals 48. 48 is an even integer and therefore the 
particular set of multiplexing factors used in the present 
invention provides one of several optimized systems. 
5 The reason this multiplexing factor product rule works well 

in a tri-level multiplexing configuration for C/A code is that 
there are three prime factors in 1023. That is, 1023 is the 
product of three prime numbers, 31, 11 and 3. Each of the three 
multiplexing factors is evenly divisible by one of these prime 
10 numbers. For example, 186 is divisible by 31 six times, 12 is 
divisible by 3 four times and 22 is divisible by 11 twice. 

Using each prime factor of the number of bits in the sampled 
/£- bit rate in one of the multiplexing factors yields two or more 
'~M different families of multiplexing configurations for C/A code 

spread spectrum receivers. In the first family, if 11 channels 
ny are desired, then' either the code repetition multiplexing factor 
:5 r or the channel multiplexing factor would have to be divisible by 
31. Although it may be desirable in certain applications to use 
\f_ 31 or 62 different code delays, there is a substantial advantage 

in making the code repetition multiplexing factor as large as 
I** possible. This reduces the number of bits required to be saved 
\Z an d processed in each segment. By selecting the code repetition 
multiplexing factor to be a multiple of 31, the number of delays 
actually used can be more easily controlled because the code 
25 delay multiplexing factor could be any multiple of 3. 

In the other convenient family, 6, 9, 12, 15 or 18 satellite 
channels are desired so that the channel multiplexing factor is 
an integral multiple of 3. This permits the code delay 
multiplexing factor to be a factor of 11 while the code 
30 repetition multiplexing factor is a factor of 31. The particular 
embodiment described in the specification above is in this 
family. 

Another constraint on the selection of multiplexing factors 
is the speed of operation of the lowest level of multiplexing. 
35 In the embodiment disclosed, the third level of multiplexing 

operates at 48f 0 . The clock speed of the hardware implementation 
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must be sufficient to permit operation at this speed. As faster 
and faster on chip components are developed, higher clock speeds 
may be used to accomplish the highest speed processing and larger 
multiples may be used. For example, with components in the high 
5 speed processing sections such as Correlator Block 110 capable of 
operation at higher rates at multiples of f 0 , such as at 96f 0 , 
the code repetition multiplexing factor could be doubled to 
produce 24 channels with 20 delays or taps or 12 channels with 40 
delays or taps or 11 channels with 6 bits and 22 taps. 
10 The system configuration may also be viewed from the 

standpoint of a time or speed magnification. Operation at the 
third multiplexing level at 48f 0 is 24 times faster than the chip 
.'S rate of the 2fo sample being processed. This amplification 

factor of 24 permits a hardware multiplexing or gate compression 
iM factor of 24. The number of gates on ASIC 102, or other devices 
ny for implementation the present invention, is reduced essentially 
: r in direct proportion to the magnification factor. All other 

factors being equal, the surface area of a chip operated at 48f 0 
if= is on the order of l/24th of the surface area that would be 
2jg required to operate at 2f 0 . Similarly, an increase in the 
\*± magnification factor to 96 would permit a reduction in the 
Jjf required chip surface real estate required on the order of almost 
half. 

The particular embodiment of the multiple level multiplexing 
25 spread spectrum receiver of the present invention which has been 
disclosed above is a GPS receiver. The same invention can be 
used for other spread spectrum signals such as wireless telephone 
signals with due consideration for the selections of multiplexing 
factors based on the bit rate of the spread spectrum code used 
30 and the environmental factors applicable to that application. 

The environmental factors for the present configuration, such as 
the pragmatic constraints on the number of channels and code 
phases, have been described above. 

Referring now to Fig. 12, a block diagram overview of a GPS 
35 receiver system 200 including a preferred embodiment of the 

digital signal processing chip 102 described above, ASIC GSP1 
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202, and a radio frequency chip, GRF1 204, combined with other 
components to form a complete receiver system according to the 
present invention . 

Associated with ASIC GSP1 202 are SRAM 206, ROM 208 and CPU 
101, interconnected by data and address busses 210 and 212 to 
provide the functions of RAM Rl 103, RAM R2 105 and Sat_Mem 186 
and other required functions described above with regard, for 
example, to Fig. 5. 

GRF1 204 is including within RF processing subsystem 214 
which received satellite signals from GPS antenna 28 and provides 
sample or GPS data 100 to ASIC GSP1 202 which returns an 
automatic RF gain control signal, AGC 216, back to GRF1 204. 
Associated with GRF1 204 in RF processing subsystem 214 are RF 
filter 218 which applies the signals from GPS antenna 28 to low 
noise amplifier LNA 220 to output of which is applied to GRF1 
204. In addition, GRF1 204 uses an outboard filter, IF FILTER 
222 as well as crystal 224. It is important to note that IF 
FILTER 222 may be a low cost, external 2-pole LC (inductance- 
capacitance) type intermediate or IF filter, rather than a more 
expensive and complex 5 or 6 pole filter for the following 
reasons. GPS receiver system 200 uses a relatively wide IF band 
followed by a decimator or digital filter, Digital Filter 118, as 
shown for example in CACAPT 104 in Fig. 5. 

In particular, the output of LNA 220 is processed by GRF1 
204 using IF FILTER 222 to produce GPS data 100 which is applied 
to CACAPT 104 in ASIC GSP1 202. Within ASIC GSP1 202, GPS data 
100 is separated into in phase and quadrature phase I and Q 
signals in I/Q splitter 106. The I signals are then applied to 
Digital Filter 118 and the Q signals are processed in the same 
manner as shown in Fig. 5 and described above. 

Referring now to Fig. 13, a series of graphs of relative 
correlation amplitude as a function of time offset are shown for 
direct path and two examples of multipath interference. These 
graphs are aligned at a time offset of zero, that is, at the time 
of arrival of the direct path signal. 

Direct path correlation function 226, in the center of the 



figure, is the idealized result of correlating a satellite signal 
received along a direct, line of sight path in the absence of 
multipath signal (s), with a replica of the C/A code modulation 
then present on the direct path signal. The peak 230 of direct 
5 path correlation function 226 is .shown at the origin to represent 
the actual time of arrival or zero code phase. In practice, this 
point may be somewhat offset due to filtering and other biases. 
Peak 230 will be taken for the purpose of this explanation as the 
punctual code phase, that is, the time of arrival of the PN Code 
10 group from a particular satellite. 

Direct path correlation function 226 may be produced for 
example by operation of Correlator Block 110 correlating the 
ji? Doppler shifted satellite signals from Doppler Block 108 by the 
output of Coder Block 112 as shown in Fig. 9 while changing the 
iM estimated code phase. In particular, direct path correlation 
jtj function 226 shows the shape of the correlation function that 
: f: would result from adjusting the code phase, in. the absence of 
" multipath interference, over a range from about one C/A code chip 
! a = width early, by a delay or time offset of -1 chip, to about one 
2^ chip width late, that is, a delay or time offset of about + 1 
j*& chip. 

\¥ The triangular -shape of direct path correlation function 226 

is, conventionally understood to result from the following 
circumstances. There will be almost no correlation between the 

25 signal received and the internally generated code when the code 
phase offset is greater than about 1 to 1.5 chips in either 
direction. As the time offset is reduced from about 1 chip to 
about zero, in either direction, the correlation increases to a 
maximum at zero offset. That is, when the code phase of the 

30 internally generated code is exactly equal (less biases, offsets 
and the effects of filtering) to the code phase of the signal as 
received, the. correlation peaks. 

A delay-locked loop is conventionally used for tracking the 
expected position of peak 230 by using a pair of early and late 

35 correlators with a fixed offset or time delay there between to 
perform early and late correlations in order to surround or 
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straddle peak 230. 

As shown in Fig. 14, residual code phase errors resulting 
from multipath interference can be detected, determined and/or 
corrected in accordance with the present invention. In 
5 particular, satellite signals are received by GPS antenna 28 and 
processed by various components as described above, as well as by 
Band Pass Filter 232, before being correlated with a code replica 
produced by PN Code Generator 234. The time offset of the PN 
code produced by PN Code Generator 234 is controlled by the delay 
10 or offset of adjustable delay 236 driven by system clock 238. 

The output of PN Code Generator 234, as offset, is applied 
to early correlator 240 for correlation with the satellite 
signals as processed by Band Pass Filter 232. The output of PN 
: <j Code Generator 234 is applied through a pair of H chip delays 242 
llf and 244 to late correlator 246, the other input of each of which 
|5] is also provided by the output of Band Pass Filter 232. As a 
: E result, the satellite signals are correlated at two points with a 
i=? fixed 1 chip delay, or separation, there between. The 
\*± correlation functions are applied to detectors 24 8 which 
2¥ evaluates a characteristic of the correlation function, such as 
the power. It should be noted that other values or 
characteristics of the correlation function, typically a complex 
number including in- phase and quadrature phase components, can 
be used in place of power measurements including amplitude 
25 measurements. 

In accordance with conventional techniques, a delay-locked 
loop is used for code tracking by adjusting the time offset of 
adjustable delay 236 so that the amplitude or power of the early 
and late correlation functions are maintained in a fixed 
30 relationship. In a preferred embodiment, the power of the early 
and late correlation functions are maintained equal by Code Phase 
Error System 250 which adjusts the code phase time offset to 
maintain this relationship. Code tracking is then performed in 
that the actual time of arrival of the code from the satellites 
35 is known to be within the one chip separation between the early 
and late correlations while their powers remain equal. 
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Referring again also to Fig. 13, when the delay-locked loop 
of Fig. 14 is shown to be properly tracking the code phase so 
that the early and late correlation amplitudes are equal, the 
relative magnitudes of early and late correlations 252 and 254 
are half of the magnitude of peak 230. That is, when the time 
phase offset is adjusted so that correlation amplitudes of equal 
value are tracked, these values symmetrically surround, in time, 
the actual time of arrival of the signals shown in the figure as 
prompt correlation 256. In other words, for a direct path 
signal, prompt correlation 256 is caused to occur midway between 
early and late correlations 252 and 254 so that prompt 
correlation 256 occurs at zero time offset, i.e. at the actual 
time of arrival of the code. As shown in Fig. 13, the amplitude 
of prompt correlation 256 is a relative amplitude assigned an 
arbitrary value of 1.0. The amplitudes of early and late 
correlations 252 and 254 have equal values of 0.5. 

As shown in Fig. 14, to cause the prompt correlation to 
occur midway between the early and late correlations, the 1 chip 
delay between early and late correlation is provided by a pair of 
H chip" delays 242 and 244.' The output of h chip delay 242 is 
applied to prompt correlator^l^e- to produce prompt correlation 
256, one half chip offset from early correlator 240, for 
evaluation by detector 248. The input of H. chip delay 244 is 
provided by H chip delay 242 so the output of H chip delay 244, 
applied to late correlator 246, is separated from the input to 
early correlator 240, by one full chip offset. The outputs of 
detectors 248 are then applied to complete the delay-locked loop. 

Multipath distortion, if present, causes the prompt 
correlation to be offset from the actual time of arrival of the 
satellite signals by an error, described herein as the code phase 
residual error. The sign of the error between the prompt 
correlation and the actual time of arrival, either leading or 
lagging, has been determined to depend upon the relationship 
between the carrier phases of the direct and multipath signals. 
When the difference in phase between the carrier phases of the 
direct and multipath signals approaches 0°, as shown for example 
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as lagging multipath correlation function 258 in the upper 
portion of the figure, the direct and multipath signals tend to 
reinforce, increasing the relative amplitude of the correlation 
products. When the difference in phase between the carrier 
phases of the direct and multipath signals approaches 180°, as 
shown for example as leading multipath correlation function 260 
in the lower portion of the figure, the direct and multipath 
signals tend to cancel, decreasing the relative amplitude of the 
correlation products . 

More importantly, the relationship between the location of 
the actual time of arrival, and the points of equal correlation 
amplitude, also changes. Without multipath, as discussed above, 
the points of equal magnitude of early and late correlations 
separated by a fixed delay are symmetrical about the correlation 
peak, that is, the actual time of arrival, so that the point 
midway there between tracked as the punctual correlation is in 
fact the actual time of arrival of the code. 

In accordance with the present invention, however, it has 
been determined that multipath interference, by reinforcement or 
cancellation, causes the points of equal amplitude early and late 
correlations to no longer be symmetrical about the correlation 
peak. For example, as can easily be seen by inspection of 
lagging multipath correlation function 258, the points of equal 
amplitude early and late correlation 252 and 254 are shifted to 
the right, that is to a positive or lagging delay, with respect 
to the points of equal amplitude early and late correlation for 
direct path correlation function 226. 

When the midway point in time offset between the early and 
late correlation is tracked, for lagging multipath correlation 
function 258, lagging prompt correlation 262 is offset in time 
from direct path prompt correlation 256 by multipath 
reinforcement interference lag error 264. That is, lagging 
prompt correlation 262 is offset from the actual time of arrival 
of the direct path signal by a positive or lagging delay time. 
Similarly, when the midway point in time offset between the early 
and late correlation is tracked, for leading multipath 
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correlation function 260, leading prompt correlation 266 is 
offset in time from direct path prompt correlation 256 by 
multipath cancellation interference lead error 268. That is, 
multipath cancellation interference lead error 268 is offset from 
the actual time of arrival of the direct path signal by a 
negative or leading lagging delay time. 

In addition, the relationship between the early, prompt and 
late correlation product amplitudes are changed by multipath 
interference. As can be seen .by inspection of lagging multipath 
correlation function 258, when the midway point in time. offset 
between the early and late correlation is tracked, for lagging 
multipath correlation function 258, lagging prompt correlation 
262 is greater in amplitude than direct path prompt correlation 
256, that is, greater than 1.0. The amplitudes of early and late 
correlations 252 and 254 for lagging multipath correlation 
function 258 are also greater than for direct path correlation 
function 226, that is, they are greater than 0.5. 

In particular, lagging prompt amplitude 270 is greater than 
1.0 and equal early and late lagging correlation amplitudes 272 
are greater than 0.5. However, as can be seen by inspection and 
as demonstrated by simulation, equal early and late lagging 
correlation amplitude 272 is greater than one half of lagging 
prompt amplitude 270. Similarly, leading prompt amplitude 270 is 
less than 1.0 and equal early and late leading correlation 
amplitudes 276 are less than 0.5. Further, equal early and late 
leading correlation amplitude 276 is less than one half of 
leading prompt amplitude 274. 

In accordance with the present invention, these 
relationships are used to determine the sign and magnitude of the 
offset errors multipath reinforcement interference lag error 2 6*4 
and multipath cancellation interference lead error 268. Code 
Phase Error System 250, shown in Fig. 14, receives as inputs, the 
correlation amplitudes (or other characteristics as determined by 
detectors 248) of the correlation products from early correlator 
240, prompt correlator 243 and late correlator 246. 

If Code Phase Error System 250 determines that the amplitude 
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of the prompt correlation performed midway between the early and 
late correlations is less than twice the amplitude of the equal 
early and late correlations, then multipath reinforcement 
interference lag error 264 exists. If Code Phase Error System 
250 determines that the amplitude of the prompt correlation 
performed midway between the early and late correlations is more 
than twice the amplitude of the early and late correlations, then 
multipath cancellation interference lead error 268 exists. 

If, however, Code Phase Error System 250 determines that the 
amplitude of the prompt correlation performed midway between the 
early and late correlations is equal to twice the amplitude of 
the early and late correlations, then no multipath interference 
error exists. 

That is, the existence of a multipath interference error may 
be detected, and if detected the sign of the error may be 
determined by comparison of the ratio of the amplitude of the 
prompt correlation to the equal amplitudes of the early and late 
correlations offset "symmetrically from the prompt correlation. 

The relative magnitude of the multipath interference error 
may be estimated in several different manners. Depending upon 
the relative amplitude of the multipath signal to the direct path 
signal, and the differences in carrier phase there between, an 
appropriate, empirically determined scale factor, such as -0.5, 
multiplied by either the sum of the early and late correlation 
amplitudes divided by the amplitude of the punctual correlation, 
or by the square root of the sum of the squares of the early and 
late correlation amplitudes divided by the square of the 
amplitude of the punctual correlation, will provide a suitable 
correction factor under most circumstances. 

In other words, a computation correction to the pseudorange 
can be made proportional to the amplitude of the (Early + Late) + 
Punctual correlations to reduce or eliminate the effect of 
multipath errors when the multipath delay is less than about 1.5 
PRN chips. 

As shown in Fig. 14, there are three uses for residual 
multipath code phase error 278, which includes both the sign and 
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estimated magnitude of the error. This error may simply be used 
in the rest of the receiver, shown as receiver processor 310, to 
computationally refine the pseudorange and therefore the position 
determination without changing the operation of the delay-locked 
loop used for tracking the code phase. 

Alternately, or in addition thereto, residual multipath code 
phase error 278 may be applied to adjustable delay 236 which 
changes the time offset of PN Code Generator 234 to control the 
offset of early correlator 240. The two H chip delays 242 and 
244 maintain the separation from early correlator 240 to late 
correlator 246 at one full chip width with prompt correlator 243 
centered there between. In this manner, prompt correlator 243 
may be made to more accurately track the time of arrival of the 
direct path signal. In addition, separation control signal 280, 
produced for example by Code Phase Error System 250, may be used 
to narrow or otherwise control the separation of the early and' 
late correlations as well as the symmetry around the prompt 
correlation to better track the actual time of arrival of the 
code. 

Further, residual multipath code phase error 278 can be used 
in multipath model 282 to enhance or provide a synthesized model 
of the interfering multipath signal (s) used, for example, for 
multipath cancellation. Replica 284 produced by multipath model 
282 may be applied as a measurement input to error correcting 
feedback loop 286 which receives the signals from Band Pass 
Filter 232 as set point input 287 to produce error signal 288 
applied to multipath model 282. Error signal 288 is used to 
control, replica 284 to reduce any differences between the replica 
and the signal received until the replica is an accurate 
representation of the multipath signals. Multipath model 282 may 
then provide additional code phase correction 290, added by 
slimmer 292 to residual multipath code phase error 278, for 
adjustment of PN Code Generator 234. 

As noted above with respect to separation control signal 
280, it may be desirable under certain circumstances to reduce 
the separation, or time offset, between early correlator 240 and 
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late correlator 246 for improved tracking. This technique of 
reducing the separation has also been used in attempts to 
minimize residual error from multipath interference by straddling 
the peak of .the correlation function. It should be noted from 
5 inspection of Fig. . 13 that the peak of the correlation function 
is easiest to select by straddling direct path correlation 
function 226 because a relatively steep slope occurs on both 
sides of peak 230, making the peak easier to distinguish. 

As the separation between the early and late correlations is 

10 reduced in the presence of multipath, however, at least one of 

the sides of the peak becomes less steep. For example, the slope 
of the lagging edge of lagging multipath correlation function 258 

:jg is substantially less steep than the leading edge thereof. 

^ Similarly, the slope of the leading edge of leading multipath 

l|l correlation function 260 is also less steep than the slope of its 

\U lagging edge. As the slope becomes less steep in the 

c 

^ neighborhood of the peak and the prompt correlation, it becomes 

more difficult to detect or straddle the peak, particularly in 
:j; the presence of noise. In accordance with one embodiment of the 
2;S present invention, the sign and magnitude of the error may first 
i a = be determined in order to reduce multipath effects before 
:5 reducing the separation for tracking purposes. 

It should also be noted from Fig. 13 that two correlation 
peaks are produced by the cancellation interference. Although it 
25 might be possible to track the wrong peak by accident, it is only 
necessary to determine if another, earlier and perhaps larger 
correlation peak exists within a small separation to cause the 
delay-locked loop to track the correct peak. 

In order to provide an accurate assessment of the multipath 
30 error, the baseline due to noise must be removed. Referring now 
to Fig. 3, the present invention permits the convenient and 
accurate assessment of the baseline due to noise as a result of 
the use of a relatively large number of correlators, such as 
Exclusive NOR gate correlators 74. In particular, a correlator 
35 74 at a substantial earlier correlation time or delay than the 
prompt tap can be used to assess the baseline due to noise. In 
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this way, the noise can be determined without interference by any 
correlation with the desired signal. 

Having now described the invention in accordance with the 
requirements of the patent statutes, those skilled in this art 
will understand how to make changes and modifications in the 
present invention to meet their specific requirements or 
conditions. Such changes and modifications may be made without 
departing from the scope and spirit of the invention as set forth 
in the following claims. 
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